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We study an inhomogeneous three-orbital Hubbard model for the Cu-substituted iron pnictides 
using an extended real-space Green’s function method combined with density functional calcula¬ 
tions. We find that the onsite interactions of the Cu ions are the principal determinant of whether 
an electron dopant or a hole dopant is caused by the Cu substitution. It is found that the Cu sub¬ 
stitution could lead to a hole doping when its onsite interactions are smaller than a critical value, 
as opposed to an electron doping when the interactions of Cu ions are larger than the critical value, 
which may explain why the effects of Cu substitution on the carrier density are entirely different in 
NaFei-xCuj;As and Ba(Fei_a;Cua;) 2 As 2 . We also find that the effect of a doping-induced disorder is 
considerable in the Cu-substituted iron pnictides, and its cooperative effect with electron correlations 
contributes to the orbital-selective insulating phases in NaFei_ 2 ;Cua;As and Ba(Fei_a;Cua;) 2 As 2 . 


I. INTRODUCTION 

Iron-based superconductors^^— have attracted great at¬ 
tention in recent years due to their considerable high 
superconducting transition temperatures and their rich 
phase diagrams, which are qualitatively similar to that 
of cuprate superconductors—!^ In cuprate superconduc¬ 
tors, substituting Cu with other transition metals is 
an effective way to gain a better insight into the ori¬ 
gin of the high-temperature superconductivity— Like¬ 
wise, many experiments were conducted to reveal the 
effects of substituting Fe with other transition met¬ 
als in some families of iron-based superconductors— “— 
Among various transition-metal substitutions, Cu dop¬ 
ing is highly disruptive to the electron structures of the 
FeAs sheets. There is lively debate on whether electron 
doping or hole doping occurs in Cu-substituted Fe-based 
superconductors— 

On the one hand, the X-ray photoelectron 
spectrapiiii^ show that Cu has a closed shell 

in Ba(Fei_a;Cu2;)2As2 and SrFe2-a;Cua;As2 compounds, 
suggesting that Cu substitution behaves like hole doping. 
On the other hand, in the 122 phase, both Cu and Ni sub¬ 
stitutions introduced a significant electron doping effect 
in Bao,6Ko.4Fe2As2 and suppressed the superconducting 
transition temperature considerably—^ Nevertheless, the 
angle resolved photoemission spectroscopy (ARPES) 
experiment on Ba(Fei_a;Cua;)2As2 predicted that a 
part of the electrons in Cu substitution preferentially 
occupied the Cu 3d states and did not behave like a 
mobile carrier— Moreover, the measurements on a 
series of NaFei_a:Cua:As single crystals showed that Cu 
substitution serves as an effective electron dopant with 
strong impurity potentialJ^ Such an observation was 
also supported by the ARPES study, which exhibited 
the increasing of the Eermi level in NaEei_a;Cua;As, 
indicating a fraction of electron doping introduced by 
the Cu dopant— 2 Up to date, the dispute on the electron 


or hole doping in the 111 and 122 phases is far from 
resolved. Therefore, further investigations are needed 
to explain the reason why the Cu substitution plays 
different roles in NaFei_a;Cua;As and Ba(Fei_a;Cua;)2As2 
compounds. 

One also notices that the random distribution of Cu 
over Fe has been confirmed by the X-ray absorption near¬ 
edge structure in the Fe K edge in Fei_a;Cua;Sei_5-^ The 
first-principles calculation also highlighted the necessity 
of including disorder effects on the iron-based supercon¬ 
ductors with transition-metal substitution^^ Therefore, 
the quantum phase transition introduced by substitut¬ 
ing Fe with Cu is complicated because Cu doping not 
only changes the carrier density but also leads to strong 
disorder effect. Additionally, the first-principles calcu¬ 
lation demonstrated that the effective on-site Coulomb 
interactions in some families of Fe-based superconduc¬ 
tors were considerably strong leading the parent com¬ 
pounds to be very close to the Mott insulating phased 
The role of electronic correlation in the renormalizing of 
electronic bandwidths^^E magnetic moments^^, orbital- 
selective Mott phase^i2^ and orbital order— in iron- 
pnictides is also recognized. More interestingly, heavy Cu 
doping drives the occurrence of a metal-insulator transi¬ 
tion (MIT)^ which rarely happens in iron-pnictide com¬ 
pounds. This naturally raises a question as to whether 
the substitution of Fe with Cu drives the system to tran¬ 
sit from a metal to an insulator, as was observed by Wang 
et al.,— and what the nature of the MIT is. 

To understand the mysterious electron or hole dop¬ 
ing and the underlying physics of the MIT tuned by 
substituting Fe with Cu in an Fe-based superconductor, 
we proposed an inhomogeneous three-orbital Hubbard 
model after reconstructing the Fermi surfaces of NaFeAs 
and NaCuAs based on the first-principles calculations. 
Here we develop an extended real-space Creen’s function 
method, which is an adequate approach for the study of 
the cooperative effect of the multi-orbital electronic cor- 
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relations and doping-induced disorder. We find that the 
increasing of the onsite interactions u of the Cu ions can 
transfer electrons of Cu ions to Fe ions, and the criti¬ 
cal point of the inversion from hole doping to electron 
doping is about Uc=5.2 eV when the interactions of Fe 
ions are 17=1.5 eV and J = U/8. It is found that the 
Cu substitutions lead to a hole doping when u <5.2 eV, 
as opposed to an electron doping when u >5.2 eV, where 
the carrier occupancy of the Cu sites drops to zero. With 
the increasing of Cu substitutions, an enhancement of the 
Anderson localization in and dyz orbitals is observed, 
whereas a Mott gap is found in the dxy orbital, indicat¬ 
ing that the cooperative effect of multiorbital correla¬ 
tions and doping-induced disorder can lead to an orbital- 
selective insulating phase in the Cu-substituted Fe-based 
superconductor. 

This paper is organized as follows. In section |nl 
we first present the electronic structures of compounds 
NaFei_a;Cua;As obtained by the first-principles calcula¬ 
tions; we then propose an inhomogeneous three-orbital 
Hubbard model to study the cooperative effects of the 
multi-orbital correlations and doping-induced disorder in 
Cu-substituted iron pnictides. The extended real-space 
Green’s function approach employed in this paper is in¬ 
troduced in section EH In section El we study the ef¬ 
fects of Cu substitution on the carrier density and dis¬ 
cuss the orbital-selective insulating phase introduced by 
the cooperative effects of multi-orbital interactions and 
doping-induced disorder. The principal findings of this 
paper are summarized in section 0 


II. LDA BANDS AND MODEL HAMILTONIAN 

We first perform the local density functional calcula¬ 
tion to understand the primary electronic structures of 
Cu-doped compounds NaFeo.gCuo.iAs, NaFeo.sCuo.gAs, 
NaFeo. 7 Cuo. 3 As, and NaCuAs, as well as the parent com¬ 
pound NaFeAs. The band structures are calculated by 
the full potential local orbital (FPL09) code^l using the 
Perdew-Wang 92 version^S of the local-density approxi¬ 
mation (LDA) for the exchange and correlation poten¬ 
tials. The orbitals of sodium (2s, 3s, 4s, 2p, 3p, 4p, and 
3d), iron (3s, 4s, 5s, 3p, 4p, 5p, 3d, and 4d), and ar¬ 
senic (3s, 4s, 5s, 3p, Ap, 5p, 3d, and 4d) ions are treated 
as valence states, while the other lower-lying orbitals are 
considered as core states. Utilizing the lattice constants 
and atomic positions obtained by experimentsthe self- 
consistency process counts a total of I04I6 irreducible 
/c-points for the fc-space integrations. 

We plot the band structures of NaFeAs and NaCuAs 
in the upper and lower panels of Fig. E respectively. The 
heavy stripes represent the weights of Fe and Cu 3d or¬ 
bitals contributed to the bands. The weights are normal¬ 
ized, that is, the summation of the weights of all valence 
orbitals at a given point of the band structure should 
be unity. From Fig. E "we find that the principal bands 
that determine the Fermi surfaces of NaFeAs comprise 


NaFeAs 




FIG. 1: (Color online) Band structures of NaFeAs (upper) 
and NaCuAs (lower) obtained by the GCA approach. Thick 
lines project out the Cu 3d orbitals in the lower figure. 


three strongly hybridized Fe 3d orbitals dxy, dxz, and 
dyz- Moreover, it is shown in Fig. Ethat the Cu 3d or¬ 
bitals primarily distribute from -2 eV to -4 eV below the 
Fermi level, suggesting that 3d orbitals of doped Cu ions 
sink below Fermi level. In addition, we find that the 
weight of those bands around the Fermi level comes from 
As Ap orbitals. The Cu 4s and Na 3s orbitals are higher 
in energy and irrelevant. We also plot the total density 
of states (DOS) of NaFeAs in Fig. E It is obvious that 
the DOS of Fe 3d orbitals at the Fermi energy is domi¬ 
nant. While, the DOS of Na 3s and As Ap orbitals near 
the Fermi level are both very small, indicating that these 
orbitals could be ignored. 

The coherent potential approximation (CPA}^!^^ is 
employed to account for the effect of substitutional 
disorder introduced by the Cu doping at Fe sites in 
NaFei_a;Cua;As. As implemented in the FPL05 code,— 
the CPA approach treats approximately the disordered 
crystal as an effective medium— In the CPA calcula¬ 
tion of the variation of electronic structure caused by Cu 
doping, the lattice parameters of NaFei_a;Cua;As (a;= 0 . 1 , 
0.2, and 0.3) are fixed at the experimental values of the 
parent compound NaFeAs, and 1496 irreducible fc-points 
are counted. 

In Fig.E we show the DOS of NaFei_a;Cua;As at a;=0.I 
and 0.3 near the Fermi level, respectively. With the in¬ 
creasing of Cu-doping concentration, the energy distribu¬ 
tion of Cu 3d orbitals does not vary significantly, clearly 
indicating that the Cu 3d orbitals always lie far below 
F,p. We do not observe MIT over a wide region of Cu- 






























3 



FIG. 2: (Color online) Density of states of NaFeAs obtained 
by the GGA approach. Three strongly hybridized Fe 3d 
orbitals d^y, d^z, and dyz determine the Fermi surfaces of 
NaFeAs. 

doping concentration, suggesting that the filling-factor 
controlled MIT could be ruled out, and thus the elec¬ 
tronic correlation or disorder driven MIT is the most rel¬ 
evant candidate. 

The band structure and orbital composition obtained 
by the first-principles calculation for undoped NaFeAs 
can be reproduced by a three-orbital tight-binding (TB) 
model. The TB model parameters are shown in Ta¬ 
ble E which are very similar to the model parame¬ 
ters of Ref. [3^]. Based on our LDA calculations of 
NaFei_a;Cua:As, we introduce an inhomogeneous three- 
orbital Hubbard model for the Cu-substituted Fe-based 
compounds, where a certain percentage of Fe sites re¬ 
placed by Cu ions are chosen at random. Therefore, the 
whole Hamiltonian consists of three parts, 

iJ = Hpe + Hc„-f Ffn.b, (1) 

where Hpe, and Hnyb represent the Hamiltonians of 
Fe ions, Cu-doping constitutions, and the hybridizations 
between Fe and Cu sites, respectively. 

The partial Hamiltonian for the electrons in dxy, 
dxz, and dyz orbitals of Fe ions can be expressed as 

a-Pcr 

i^Fe,cit<T i^Fe,CK 

+ E E ({7 JS(ja')n'iaa'nil3a' 1 (2) 

i^Fe a</3,(T(T^ 

where iciaa) creates (annihilates) an electron with 
spin projection a for orbital a. {dxz,dyz, or dxy) of an iron 



Energy (eV) 


FIG. 3: (Color online) The density of states of compounds 
NaFei_a;Cu 2 ,As at a;=0.1 (a), 1=0.2 (b), and a;=0.3 (c) within 
the coherent potential approximation. 

site i, and riiacr = cl^a^iacr represents the corresponding 
electron number operator. represent the hopping 

integrals of electrons between nearest-neighbor (NN) or 
next-nearest-neighbor (NNN) iron sites i and j within 
orbitals a and /?.— Axy = fJ,xy — fJ^xz is the energy dif¬ 
ference between the dxy and degenerate dxz/dyz orbitals. 
More details regarding the definition of kinetic energy in 
a Hamiltonian ([2]) can be found in appendix U and 
U' are the on-site intra-orbital and inter-orbital Coulomb 
interactions on iron sites, respectively, and J = Jz rep¬ 
resents the Ising-type Hund’s couplings. In our calcula¬ 
tions, we ignore spin-flip and pair-hopping processesi^^— 


TABLE I: Parameters for the hopping integrals of the three- 
orbital Hubbard model.— The energy unit is electron volts. 
The quantities are defined in appendix [Al 


Ti 

T2 

Ts Ti 

n n 

TV Tg Axy 

0.15 

0.15 

-0.12 0.06 

-0.08 0.1825 

0.08375 -0.03 0.75 


As shown by the first-principles calculations, the hy¬ 
bridization is mainly between the Fe 37 orbitals and the 
Cu 4s orbitals, but it is very weak. The Cu 37 orbitals 
distribute from -2 eV to -4 eV below the Fermi level, 
unlike the Cu ions in cuprate high-Tc superconductors. 
When an Fe ion is substituted by a Cu ion in Fe-based 
superconductors, we need to consider all three t 2 g or¬ 
bitals of a Cu ion instead of the single dx^-y^ orbital. 
Therefore, the Hamiltonian for the Cu substitutions is 

Hcx, = — E! 

7'r 

- V-yuf^^+U Y EtEt 

iGCuy’ycr i^Cu,'y 
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+ Y1 (3) 

iGCu 7 <A,crcr' 

where (di'^a-) creates (annihilates) an electron at a 
copper site i'. are the chemical potentials of different 
orbitals 7 . The intraorbital hopping integral between 
two copper sites is indicated as Uiji, and the intraorbital 
and interorbital interactions, as well as the Hund’s rule 
couplings, are represented by u, u', and j for copper sites, 
respectively. 

As our calculations show that the intraorbital interac¬ 
tions u play a principal role, only the effect of u is taken 
into account. That is, we choose u' = u and j = 0 in our 
discussion of Cu substitution effect in the iron pnictides. 
Because the Cu Sd orbitals locate about 3 eV below the 
Fermi level, the intraorbital hopping integrals tij are pre¬ 
dicted to be quite small. Owing to the restriction of the 
precision limitation of our calculations for the disordered 
cases, we approximate tij in the range of 0 to 0.03 eV for 
the NN hopping integrals between Cu ions. In our calcu¬ 
lations, we choose t=0.01 eV for the NN hopping terms 
but t = 0 for other long-distance hopping terms. We find 
that the value of t only has an extremely-limited effect 
on the bandwidth of the Cu impurity band, and its effect 


on the DOS of Fe is negligible when the concentration of 
Cu substitution is less than 0.2. 

Moreover, the hybridization of Fe and Cu ions can be 
expressed as 

Hnyb = - ^^J(4aadJ'ra + il■C.). (4) 

iGFe,OL j^Cu^'jcr 

Because of the very weak effects of the long-distance 
terms, we only need to consider the NN hopping in¬ 
tegrals d between Fe and Cu ions, which are also as¬ 
sumed as 0.01 eV. Just the same as t, the effect of t' 
on the DOS of Fe ions is also negligible for the case with 
X < 0.2. We study the inhomogeneous three-orbital Hub¬ 
bard model Eq. o by developing a real-space Green’s 
function method, which is described in the next section. 


III. REAL-SPACE GREEN’S FUNCTION 
APPROACH 

For a three-orbital Hubbard model, the real-space 
Green’s function of a square lattice with N = sites 
is expressed as a 3A x ZN matrix 


Gii 

••• Gi,h-i 

FL 

Gi^h+i 

Gim ^ 


Gh-1,1 

Fhi 

Gh+1,1 

■ ■ ■ Gfi—l^h—l 

Fh,h-i 

■ • • Gh+i,h-i 

^h-l,h 

Dhh 

p* 

Gh-l,h+l ■ • ■ 

Fh^h+i 

Gh+l,h+l ■ • • 

Gh-1,M 

FhM 

Gh+1,M 


Gmi 

■ ■ ■ GM,h-l 

p* 

^ Mh 

GM,h+l ■ • ■ 

Gmm / 

M^3N 


(5) 


where h represents the sites occupied by Gu ions. Gij, 
Fij, and Dij are 3x3 matrices with elements defined as 

Gtf = I 

Fr/ = I 

Dtf = I 4J). (6) 

Each element of the Green’s function matrices is obtained 
by the equation of motion,— for example, 

(lO -\- ^a){{Ciacr \ = SijSajS ~ ^ ^ tib{{dbjcr \ ^jpu)) 

6 GCn ,7 

~ 'y Tibamiicbma \ Cjjj^)) 

b^Fe,m 

+U{{ 

‘^iaaCiaa \ )) F U ^ '^ {{nila-Ciaa- \ cj^g.)) 

l^a 


F{U J) 'y^^iinuo-Ciaa I (7) 

l^Oi 

from which come the second-order Green’s functions, i.e. 
{{niiifC^aa I and {{nn^Ciaa I on the right 

side of Eq. 0 . Similarly, the second-order Green’s func¬ 
tions can be obtained by their equations of motion, where 
we will find the third-order Green’s functions. Having 
an analogy with the Hubbard-I approximation^ we in¬ 
troduce a proper decoupling process to make up a self- 
consistent loop for the calculations of inhomogeneous 
Green’s function G. More details can be found in Ap¬ 
pendix |B] 

We thus obtain a complete and solvable set of equa¬ 
tions for all elements of the Green’s function matrix G in 
Eq.([5D, 
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Ma{{Ciacr \ 

Mo,{{Ciaa I d^jXa)) 

(w + l'j){{dija I Cjfia-)) 

(w + V'y){{dija I d^jXa-)) 


^a^ij^a0 (-^q d^a){ ^ ^ ^ibctmU^'bma \ ^ ^ I 

b^Fe,m bGCu,'y 

U'-J 


+- 


A„ 


' y ] ( y ] Tiblm{c\j^^Ciaa) d- ^ ^ t^b{dlj^Ciaa^)){{Cila- |cj|^^)), 


l^a. bGFe,Tn b^Cu,’y 

— {Na — Ka){ 'y^ Tn,am{{cbma \ ^^jXa)) ^ibii^b^a \ d^jXa))) 


b^Fe,m 


b^Cu,’y 


U'-J 


A„ 


y ] ( y ] T'iblm.icljj^^Ciatr) + ^ ^ t.if^{dl^^^Cia(j)){{Cil(7 \ '^jAcr))> 


6 GCn ,7 


^ C((Caa. ISV))), 

^ bGCu aGFe,a 

(1 + J fe-- E I^IaJ)- E I^]aJ)), 


b^Cu 


a^Fe,a 


r 


( 8 ) 


with 

M„ = 


iV« = 

— 

Ba = 

= 


W + d'ct + 2 y ^ tif^{Ba{djj^^Ciacr) 

6 GCn ,7 

— J / jt \ \~^ / jt \ \ 

■I 4 ^ i^b'fa^il-’^/ + ~A~ \dbjgCila)), 

, \ \~^ / \ / \ 

1 H“ ^a.\^ia.a) 7 / ^ {'f^ila/ “t“ “7 / ^ y^ila} ■) 

■^ot , , -^Q; , , 

Z^a Z^a 


CJ • 


l^OL U ^ ^ ij^isa) H" J ^ ^ {j^isa ) ; 
S^Ol,(7 s^oc 


U'-J 




(9) 


l^Ct 


Therefore, the real-space Green’s function matrix ^ 
can be obtained self-consistently by the matrix relation¬ 
ship 


^ = JA, 


( 10 ) 


where ^ and jV are also 3A^ x 3A^ matrices that consist 
of the correlation functions, which can be obtained by 
the spectral theorem {AB) = —\ f(co)Im((A | B}}. 
The real-space Green’s functions can be transformed into 
the momentum-space by using the Fourier transforma¬ 
tion, G^^(a;) = where k 

denotes the wave-vectors in the unfolded Brillouin zone 
(BZ)j^ The one-particle spectral density is obtained by 

A(k, to) = -ilmGij.(a;), (11) 

where the matrix of the momentum-space Green’s func¬ 
tion Gj,(u;) is diagonalized to include the strong hy¬ 
bridization of dxy, dxz and dyz orbitals of Fe ions. In the 


next section, we study the effect of disorder introduced 
by substituting Fe with Gu in Fe-based superconductors. 

To study the substitution effect of Cu ions in a finite 
square lattice, a portion of sites chosen randomly should 
be assigned to the Cu impurities. For example, if the con¬ 
centration of the Cu substitution is x = 0.04, we choose 
randomly 16 sites for Cu ions in a IV = 20 x 20 lattice, 
and all the other sites are for Fe ions. For a certain dis¬ 
order configuration, we calculate the real-space Green’s 
function as described above to obtain physical properties, 
such as the local density of states (LDOS) at site i for 
orbital a, 


Pa(r„a;) =--ImG““(a;). (12) 

TT 

To find the averaged DOS of the whole system, we need 
to calculate the LDOS of different disorder configurations 
and then determine the sample-averaged values. The 
DOS of Cu and Fe ions can be expressed respectively 
as 




N, ^ Wp. 


ieFe 


JCu) = 1 ^ p«(D,a;), (13) 


Ns ^ Ncn 


i^Cu 


where Npe and Ncu represent the total numbers of the Fe 
and Cu ions, respectively. is the number of disorder 
conhgurations, and we consider more than one hundred 
disorder samples in the following calculations. The effect 
of disorder could be studied accurately because our nu¬ 
merical method is just the same as the real-space exact 
diagonalization method^ when all the terms of interac¬ 
tions are omitted. 
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FIG. 4: Spectral density along the high-symmetry directions 
in the unfolded Brillouin zone with different intraorbital in¬ 
teractions: t/ = 0 (a), t/=l eV (b), and C/=1.5 eV (c) for 
undoped NaFeAs, where the Fermi energy is indicated by 
a dashed line. The Hund’s rule coupling is J = U /8, and 
U' = U — 2J is satisfied. 


IV. COOPERATIVE EFFECT OF 
INTERACTIONS AND DISORDER 


The strong correlation and disorder, as described 
above, have been recognized as two possible fundamental 
origins that can drive the MIT. In the correlated fermion 
systems with the presence of disorder, the competition 
between Anderson and Mott-Hubbard MITs is still a 
theoretical challengers Iron-based superconductors are 
typically multiorbital correlated materials. The ARPES 
measurement reported strong interactions between elec¬ 
trons in NaFeAsiSS Furthermore, in NaFei_a;Cua;As, a 
strong impurity potential is introduced by substituting 
Fe with Cu randomlyrS In this section, to elucidate the 
Cu substitution effect in iron pnictides doped with cop¬ 
per, we use the newly developed real-space Green’s func¬ 
tion method described above to study the cooperative ef¬ 
fect of the multi-orbital interactions and doping-induced 
disorder. 


A. Multi-orbital Interactions in nndoped NaFeAs 

To estimate the multiorbital interaction parameters of 
Fe ions, such as the intraorbital interactions U, interor¬ 
bital interactions U', and Hund’s rule coupling J, we fit 
our numerical data of the band structures to some ex¬ 
perimental and theoretical findings, as shown in Fig. U) 
Because the undoped Fe-based superconductors are indi¬ 
cated to have a filling of approximately two-thirds based 
on the band structure calculationswe determine the 
Fermi energy by using a constraint of four electrons per 



FIG. 5: (Golor online) Fermi surfaces of the undoped three- 
orbital Hubbard model with U = 0 (a), U=1 eV (b), and 
17=1.5 eV (c). In an area of the unfolded Brillouin zone, 
there are two hole pockets around the F point (0, 0) and one 
electron pocket around the X point (tt, 0). The dashed lines 
show schematically the Fermi surfaces of the ARPES study 
on the parent compound of the 111 phase.— The model 
parameters are the same as those in Fig. 21 


Fe, 


1 

ntotai = -^ / /(a;)ImG““(w)dw = 4, (14) 

ia 

where /(w) is the Fermi-Dirac distribution function. 

Being consistent with the noninteracting spectral den¬ 
sity of Ref. [22 ], all three t 2 g bands across the Fermi sur¬ 
face with bandwidths of £>1 = 1.95 eV, £>2=0.9 eV, and 
£>3=1.5 eV, respectively, are shown in Fig. UJa) from 
high energy to low energy. A remarkable effect on the 
spectral function is found when the multi-orbital inter¬ 
actions are taken into account, especially when all three 
bands reduce their bandwidth with the increasing inter¬ 
actions. According to some previous theoretical studies 
of the iron pnictides,— Hund’s rule coupling is pre¬ 
dicted to be J = C//8 and the relation U' = U — 2J is 
also satisfied. As shown in Fig. UKc), the bandwidths of 
the three orbitals are £>i=l eV, £>2=0.7 eV, and £>3=1 
eV, respectively, when U is increased to 1.5 eV. Approx¬ 
imately, the on-site Coulomb repulsion renormalizes the 
bandwidth of the three bands across the Fermi surface 
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by factors of 2, 1.3, and 1.5, respectively, when ?7=1.5 
eV. 

Compared with the ARPES measurement rea¬ 
sonable spectral functions and Fermi surfaces are repro¬ 
duced qualitatively for undoped Fe-based superconduc¬ 
tors NaFeAs and LiFeAs by using the numerical calcula¬ 
tions of the undoped three-orbital Hubbard model. That 
is, when U= 1.5 eV, the bandwidths of the spectrum 
are about the half of the bandwidths obtained by the 
LDA, in good agreement with the ARPES experimen¬ 
tal observation!^ Thus, the choice of the intraorbital 
interactions as U=1.5 eV is reasonable in the following 
numerical calculations of the Cu substituting effects in 
NaFei_a;Cua;As. 

Apart from the narrowing of the bandwidths, a signif¬ 
icant shift of the Fermi level is also found accompany¬ 
ing the increasing interactions. The three-orbital tight- 
binding model exhibits two d^^z/dyz-dominant hole pock¬ 
ets around the P point and one da;y-dominant electron 
pocket around the X point. The change of the interac¬ 
tions can only transfer electrons from dxz/dyz orbitals to 
dxy orbital. With increasing U, a clear evidence of the 
increasing of the outer hole pocket is shown in Fig.[5l Ac¬ 
cordingly, the sizes of the electron pocket and the inner 
hole pocket also expand with the increasing intraorbital 
interactions, but are not as significant as the size of the 
outer hole pocket. Our numerical results are qualitatively 
compatible with the LDA-f DMFT study on the effects of 
interactions in LiFeAs j^^i^^ 

For the convenience of comparison, we also display 
schematically the Fermi surfaces of ARPES study on the 
parent compound LiFeAs .^^’*^^ As shown in Fig.[5l[a) and 
Fig. Eljb), both the electron and the outer hole pockets 
obtained in the numerical simulations are smaller than 
the corresponding experimental results in either case, 
i.e. when U=0 or U=l. As shown in Fig. El^c), we can 
find the best fitting of the Fermi surfaces obtained nu¬ 
merically with the experimental results, suggesting that 
the reasonable value of the onsite interactions of Fe ions 
should be 17=1.5 eV. Here we need to emphasize that, 
our theoretical model is applicable to both the 111 and 
122 phases. As indicated by the ARPES experiments^, 
in LiFeAs and BaFe 2 As 2 , the sizes of the electron and 
hole pockets are comparable. 


B. Effect of Cu substitution on carrier density 

Some experiments concerning the effect of Cu substitu¬ 
tion are mainly carried out in the 111 -class and 122 -class 
iron-based superconductors. In the 111 phase, such as 
NaFei_a;Cua;Asii and LiFei_a;Cua;As,— it is found that 
Cu substitution enlarges the electron pocket and shrinks 
the hole pockets, implying that some extra electrons are 
brought in. Meanwhile, in the 122-phase, the experi¬ 
mental results have not reached agreement yet, but the 
prediction of hole doping, introduced by Cu substitution 
in 122 compounds, is better received. The 3d states of a 
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FIG. 6: (Color online) The average density of states of iron 
(solid line) and copper (dashed line) sites over 100 disorder 
configurations for different interactions of Cu ions: u=3.6 eV 
(a), u=5.1 eV (b), and u=5.4 eV (c). The intraorbital interac¬ 
tions of Fe ions are U — 1.5 eV, and the concentration of the 
Cu substitution is x = 0.04. The lattice size is A = 20 x 20, 
and the total number of disorder samples is = 100. 


Cu ion in Ba(Fei_a;Cua;) 2 As 2 were found at the bottom 
of the valence band, forming a localized shelh ^^i^^ 
As a result, the Cu substitution should result in a hole 
doping in the 122 phasei ^^d^ 7 ^i ^6 Although the volume 
of the electron Fermi surface is smaller than the predicted 
value of the rigid-band model, effective electron doping 
introduced by the Cu substitutions was also predicted by 
other experiments ] 
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FIG. 7: (Golor online) The average carrier occupancy on each 
copper site as a function of the interactions u. Inset: The sub¬ 
stitution dependence of nc-u for different interactions: m= 5.1 
eV (a), and m= 5.4 eV (b) when t/=1.5 eV and J = t//8. 
The lattice sizes are A = 20 x 20, and one hundred disorder 
configurations (As = 100) are considered. 

Next, we discuss the distinct effects of Cu substi¬ 
tutions on the carrier density in NaFei_a;Fea;As and 
Ba(Fei_a;Cua;) 2 As 2 . A first-principles calculation has 

































indicated that the screened Coulomb interaction at Fe 
atoms of 122 families is less than that of Fe atoms of 
111 families for approximately 0.35 eV.— As a result, for 
the two types of iron-based superconductors, the onsite 
interactions of the Cu impurities should be different ac¬ 
cordingly. Because of the comparative localization of the 
Cu orbitals, the effective interactions of Cu ions should be 
considerably large. By contrast, the wave functions of Fe 
orbitals that cross the Fermi level are extended, introduc¬ 
ing small effective interactions of Fe ions. For simplicity, 
in our calculation we assume that the interactions at Cu 
sites (u) are in direct proportion to the interactions at Fe 
atoms ([/), with an approximate relation u ~ 4C/. Con¬ 
sidering that in Fe-based superconductors, the value of 
U should be 1 eV< 17 <2 eV, it is reasonable to predict 
the interactions of Cu ion as u ^ 4-7 eV. The above esti¬ 
mated parameters have been proved to be appropriate by 
a combined valence band photoemission and Auger spec¬ 
troscopy study of single crystalline Ba(Fei_a;Cua;) 2 As 2 , 
where the interactions in Fe and Cu ions are found to be 
(1.4±0.6) eV and (7.5±0.4) eV, respectively^^ Assuming 
that u ~ 417, the deviation is approximately 1.4 eV for 
the Cu interactions of the 111 and 122 phases. 

Effects of the onsite interactions of Cu ions on the DOS 
are shown in Fig. [6l where the concentration of Cu sub¬ 
stitutions is fixed at a: = 0.04. For simplicity, the interac¬ 
tions of Fe ions are also fixed as 17=1.5 eV and J = 17/8. 
When u = 3.6 eV, we find very narrow upper and lower 
Hubbard bands of Cu impurity that locate at w = 0.7 
eV and w = —2.4 eV respectively, as shown in Fig. |6l[a). 
Furthermore, the DOS of lower Hubbard band is con¬ 
siderably larger than that of the upper Hubbard band, 
suggesting that most of the electrons in Cu ions locate far 
below the Fe bands when u = 3.6 eV. We also find that 
the DOS of the copper sites shifts from the lower Hub¬ 
bard band to the upper band with increasing u as shown 
in Fig. mb) and (c). In addition, the upper Hubbard 
band moves upwards, such as the upper Hubbard band 
raises to a;=4.3 eV when w=5.4 eV, which is far above 
the Fe bands. Our numerical results demonstrate that 
electrons in the Cu band can transfer to the Fe bands 
with the enhancement of the interactions of Cu ions. 

As shown in Fig. [3 the carrier occupancy on Cu sites 
is obtained by integrating the corresponding DOS of Cu 
ions. We find that the average electron occupancy on Cu 
sites decreases with increasing u, especially the curve of 
ricu vs M, which drops sharply in the vicinity of u=5.2 
eV. From the perspective of the carrier density on the 
Fe sites, the Cu substitutions lead to a hole doping when 
u <5.2 eV because most electrons rest at the Cu orbitals. 
However, it is entirely different when u >5.2 eV, where 
the carrier occupancy of the Cu sites drops to zero, lead¬ 
ing to electron doping. 

To make it more clear, the inset figures in Fig.[7]show a 
reversal tendency of the dependence of ricu on x for differ¬ 
ent u. It is shown that the electron occupancy of each Cu 
site decreases with increasing x within the hole doping 
region with m= 5.1 eV. On the contrary, the occupancy 


increases with increasing x in the electron doping region 
with m= 5.4 eV. The ARPES experiment on NaFeAsii 
indicates that the total number of electron carriers intro¬ 
duced by Cu atoms increases when the Cu substitution 
increases from 0.019 to 0.14. However, the average num¬ 
ber of electrons in Cu atom decreases significantly from 
1.9 to 0.64 with increasing doping concentration. Erom 
the inset in Fig. [3b), we find that the average electron 
occupancy on a single Cu site also increases with increas¬ 
ing X, implying that less electrons can be pushed from Cu 
sites to the iron sites. This tendency is consistent with 
the ARPES experiment. 



FIG. 8: (Color online) The variation of E/, AE/ = Ef — Efo, 
vs the interactions u of Cu ions, where Efo is the Fermi energy 
in undoped compound. The lattice parameters are the same 
as those in Fig. 0 

In Fig. m we plot the dependence of the variation of the 
Fermi energy Ey on the interactions u. In the weak inter¬ 
action range of u, the Fermi energy of a Cu-substituted 
compound is less than that of an undoped compound, 
being consistent with the hole doping feature. However, 
when u >5.2 eV, AEy becomes positive, which means an 
electron doping appears. Our calculations demonstrate 
that the Coulomb interactions u affect the average car¬ 
rier density on iron sites and the Fermi energy of the 
whole system by pushing the electrons out of copper or¬ 
bitals. For this reason, the different effects of Cu substi¬ 
tutions on Fermi energy and carrier density in 111 and 
122 phases can be understood because the correlations of 
the 122 families are smaller than the correlations of the 
111 families. 


C. Localization effects of Cu substitutions 

Aside from the change in carrier density, another re¬ 
markable effect of the Cu substitutions is the doping- 
induced disorder. Adding the Cu substitutions, the 
bands are most strongly modified near the Fermi level. It 
is shown that both two hole pockets around the F point 
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(b) 


(c) 


00-jI 


FIG. 9: (Color online) (a) Spectral density of Cu-doped 
NaFeAs with Cu impurity concentration x = 0.04. (b) and (c) 
show the density of states of the degenerate d^z/dyz and d^y 
orbitals for x = 0 and x = 0.04, respectively. The Fermi en¬ 
ergy level is represented by the black dashed line. The lattice 
size is = 20 X 20, and the number of disorder configurations 
is Ns = 100. The other parameters are t/=1.5 eV, J = t//8, 
and u=6 eV. 


and one electron pocket around the X point are very sen¬ 
sitive to the Cu substitution. Apart from the shrinkage 
of the bandwidth, the most dramatic change in the band 
structures in NaFei_a;Cua;As is that the spectral density 
near the Fermi surface is strongly suppressed, as shown in 
Fig.EKa). As a result, the compound exhibits insulating 
or bad metallic behavior. 

Efros and Shklovskii have demonstrated that the in¬ 
teractions between the localized electrons in a disordered 
system can create a Coulomb gap in the DOS near the 
Fermi levelj^ As a result, the appearance of the zero-bias 
anomaly at the Fermi energy is seen as proof of the lo¬ 
calization of electronic states in the correlated systemsj^ 
Just as expected, a soft gap can be found at the Fermi 
level in the DOS of both the degenerated dxz/dyz orbitals 
and the dxy orbital when x = 0.04, as shown in Fig. IHKb) 
and (c). The characteristic of the DOS of undoped com¬ 
pound is entirely different because there is no soft gap 
at Fermi level for all three orbitals, as shown by the red 
dashed lines in Fig. |9Kb) and (c). Moreover, the Efros- 
Shklovskii gap can also be found in Fig. [SI Therefore, 
we believe that the disorder effect in NaFei_a;Cua;As is 
considerably strong. 

The iron-based superconductors are predicted to be 
close to the MIT border but are in the metallic side4^ 
There are three possible origins for the MIT in the corre¬ 
lated system with disorder: (I) Mott-Hubbard MIT: the 
on-site Coulomb interactions split the energy band near 
the Fermi level into higher and lower Hubbard bands; 
(2) band occupation-driven MIT: changing the carrier 
density or bandwidth to make the conducting band fully 
filled; (3) Anderson MIT: electrons are localized in the 
real-space by Anderson localization effects introduced 
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FIG. 10: (Golor online) The evolution of the density of state 
of the dxy orbital (the left panel) and the dxz/dyz orbitals (the 
right panel) for different Gu impurity concentrations: x = 0 
((a) and (b)), x = 0.02 ((c) and (d)), and x = 0.08 ((e) and 
(f)) when J = U/8 (dashed line) and J = U/5 (solid line), 
respectively. The other parameters are: U=1.5 eV and u=6 
eV. The lattice sizes are A = 20 x 20, and the number of 
disorder samples is As = 100. 


by disorder. Experiments found Cu substitution-driven 
MIT in NaFei_a;Cu 2 ,As, but the mechanism of the MIT 
is still on debatable. 

When 17=1.5 eV and J=U/8, it is found that Ander¬ 
son MIT happens with the increasing of the Cu substitu¬ 
tions, owing to the effect of disorder introduced by the Cu 
dopant in Ee-based superconductors. Besides, there is no 
Mott gap for both the dxy and dxz/dyz orbitals at Fermi 
level when J = U/8, as shown in Fig. (TUI In general, the 
multi-orbital interactions, especially the Hund’s rule cou¬ 
plings are likely to cause orbital-selective transitionsj^ It 
has also been predicted in some theoretical studies that 
a wide range of J/U is possible for the iron-based su¬ 
perconductors. For example, the proper range of J/U 
couplings was predicted as 0.1< J/U <0.33, which was 
obtained by comparing the theoretical results with neu¬ 
tron scattering and ARPES experiments)^ Therefore, we 
also study the effect of the Hund’s rule coupling on the 
MIT in the Cu-substituted iron pnictides. An orbital- 
selective insulating phase is found when we increase the 
Hund’s rule coupling to J = U/b and the on-site inter¬ 
actions of Cu ions to M = 417=6 eV. By observing the 
evolution of the DOS of the three 3d orbitals, we find 
that the properties of the insulating phase of the dxz/dyz 
and dxy orbitals are different. As shown in the left panel 
of Eig. 1101 a zero-bias anomaly dip can be found with 
increasing Cu substitution in the DOS of the dxz/dyz or¬ 
bital, but no hard gap can be obtained. Therefore, the 
lower-energy states at the Fermi level are Anderson local¬ 
ized as a result. Conversely, as shown in the right panel 









































of Fig. [TOl a hard Mott gap is opened in the Fermi level, 
indicating that the d^y orbital is in the Mott insulating 
state. In addition, when J = C//5, the carrier type also 
changes from a hole doping to an electron doping with 
the increasing of the Coulomb interaction u of Cu ions. 

The effect of disorder on the conductance of Fe-based 
superconductors with Cu doping is studied by using the 
Kubo formula )^^'^^ which can be simplified, when w = 0, 
as a 3 X 3 matrix a for the three orbital Hubbard model, 

a = -l,j -1) 

TT 

TlrnGjO' - 1,/ - l)ImGj(j',j) 

-ImG,(j,/ - l)ImGi(j',j - 1) 

-ImG,(j - l,/)ImGj(/ - l,j)]. (15) 

The trace is over the sites perpendicular to the current 
direction (fc direction), (i, j) denote the coordinates of 
lattice sites in the ribbon, where the periodic boundary 
conditions should be applied only to the i direction, and 
j and j' are chosen to be on opposite sides of the sam¬ 
ple along the current direction. The conductance of the 
whole system can be obtained by adding the conductance 
of all three orbitals, which correspond to the three diago¬ 
nal elements of the conductance matrix a. As mentioned 
in Sec. imi the averaged conductance is extracted by 

S = ( 16 ) 

m 

where is the conductance obtained by Eq. (fTSl) for 
the m-th disorder configuration, and the resistance could 
be obtained directly by p = 1/d. 

To compare with the experimental results, in Fig. 1111 
we show the temperature dependence of the resistance 
for different Cu concentrations. It is obvious that the 
compound behaves as an insulator when the density of 
the Cu substitution is over x — 0.1. We find that the 
MIT in NaFei_a;Cua;As is a transition from a metal to an 
orbital-selective insulating phase, which is caused by the 
cooperative effect of multi-orbital interactions and the 
doping-induced disorder. With the increasing of Cu sub¬ 
stitutions, the effect of disorder significantly enhances, 
leading to strong Anderson localization in the dxz and 
dyz orbitals. Meanwhile, a Mott gap occurs in the dxy 
orbital. As a result, the resistance increases with increas¬ 
ing doping concentration both for the whole system and 
the individual orbitals, in good agreement with the ex¬ 
perimental results)^ 


V. CONCLUSION 

We present the study of an extended real-space Green’s 
function method combined with density function cal¬ 
culation on the cooperative effects of multiorbital cor¬ 
relations and disorder in Cu-substituted iron-based su¬ 
perconductors. The inhomogeneous three-orbital Hub- 
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FIG. 11: (Color online) The temperature dependence of the 
resistance with increasing Cu substitution from a: = 0.1 to 
X = 0.16. (a) Experimental data taken from Ref. [60]. (b) 
Numerical result of the resistance for the whole system. Addi¬ 
tionally, the resistance of each orbital for d^z/dyz (c) and d^y 
(d). Each data point is obtained by averaging for one hundred 
disorder configurations. The lattice sizes are N = 20x20, and 
one hundred disorder configurations {Ns = 100) are consid¬ 
ered. The other model parameters are the same as in Fig. 1101 


bard model has been employed to describe the substi¬ 
tution of Cu with Fe based on the electron structures 
of NaFei_a;Cua;As compounds. We find that the type 
of charge carriers introduced by the Cu substitution, 
namely, electrons or holes, is virtually determined by the 
effective interactions of the copper atoms. Due to the 
different environments surrounding the Cu substitution, 
the interactions of Cu ions in 111 phase are larger than 
that of 122 phase. As a result, a hole doping is introduced 
by the Cu substitution in Ba(Fei_a;Cua;) 2 As 2 , as opposed 
to an electron doping in NaFei_a;Cua;As. Moreover, an 
orbital-selective insulating phase is observed where the 
dxz and dyz orbitals are Anderson localized, while the 
dxy orbital is in the Mott insulating phase. These results 
explain the experimental observations in NaFei_a;Cua;As 
and Ba(Fei_a;Cua;) 2 As 2 . 
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TABLE II: The NN and NNN hopping integrals of the electrons in the three orbitals d^y, d^z and dyz of Fe ions. (—1)'*' 
indicates that the parameters change sign along the site locations. 
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Te 

Tr 
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T7 
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(yz,xz) 
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-Ta 

Ta 

-Ta 

Ta 


Appendix A: Hopping integrals in the three-orbital 
tight-binding model 

As shown in Figs. mm the local density function calcu¬ 
lation predicted that the Fermi surface of NaFei-^Cu^^As 
composed only three Fe 3d orbitals {dxy, dxz and dyz), 
whereas Cu 3d orbitals contributed from -4 eV to -2 eV 
under the Fermi level. (Eq. (2). Daghofer et al. have 
constructed a three-orbital Hamiltonian for the pnictides 
following the Slater-Koster procedurej ^^i^^ In keeping 
with their model, we mainly consider the NN and next- 
nearest-neighbor (NNN) hopping integrals of electrons in 
the three orbitals dxy, dxz, and dyz of iron ions. The def¬ 
initions of the different types of hopping terms Tiajp in 
the tight-binding model are displayed in Table HIl 

The hopping integrals between the dxy and dxz/dyz 
orbitals contain factors (—1)* that arise from the two- 
iron unit cell of the original FeAs planesAs a result, 
those interorbital hopping terms change sign for NN iron 
sites. Moreover, the NN intraorbital hoppings between 
dxy and dyz along the x direction, and between dxy and 
dxz along the y direction, as well as the NN interorbital 
hopping between dxz and dyz orbitals, are all omitted. 

Apart from the hopping terms, we need to add a term 
in the Hamiltonian for the energy splitting between the 
dxy and degenerate dxz/dyz orbitals, l\xy, which is cru¬ 
cial to providing the proper Fermi pockets when fitting 
with the local density approximation results. 


Appendix B: Decoupling approach 

All approximations that we used to decouple the third- 
order Green’s functions to obtain Eq. m are listed as 


follows: 



{{fT'iaa' Cbma 


^ (niacr')((cbmcr | 

( dtja 

1 - 

{niaa') {{db'fcr \ C^p^)), 


1 - 

^ (^bmd-^io,S-}((Ciacr | c]]^^)). 


1 " 

{^iaa^bras) {{Ciaa \ c]]^^)). 



1 cl 

1 S/3o 

)) ^ 

{^bma^'i-ls) {{^iaa I 

b'jpfj))’’ 


1 cl 

1 S/3o 

)) ^ 

^ {(^ila^bma) {{Ciaa \ 



1 cl 

1 ^jl3a 

)) ^ 






— (^bma- ((Ciid 

1 sV))’ 

( {Cj^ifjCbrrKT ^ioLO 

1 cl 

1 S 730 

)) ^ 

{^ila^bmu) {{Ciaa | 





— ( {^bmcr 

1 sV))’ 

^icta 

1 cl 

1 

)) ^ 

^ {{Ciaa 1 

''j/3d')). 


1 cl 

1 S/3o 

)) ^ 

{^iaS^b-fo-) {{Ciaa | 

b^jpcr))^ 

^^db^^C-ila C-icta 

1 cl 

1 S/3o 

)) ^ 

^ {{c-ian \ 



1 cl 

1 

)) ^ 

^ i^lla^b-fa) {{Ciatj \ 

S'/3d)). 

{{df^^^Cila-Cicxa 

1 cl 

1 

)) ^ 

^ ('^b-yo-Czid) ((Ciad- 1 





— (b^bycT^^aa) {{Cila 

1 b^jPa))^ 


1 cl 

1 S/3o 

)) ^ 

^ (b^jlcrb^b'/a) ((Cia/7 | 

c]/3d)) 




— (b'ila-b'iaa) ( {db-fa 

1 b^jPa))^ 

'^ioiaC-icta 

1 cl 

1 S/3o 

)) ^ 

^ {{‘^iaa’Ciaa 

c]/3d)). 


1 cl 

1 S 730 

)) ^ 

^ ^icx<7 

1 





where a' and 

represent 

a or a. 



In the above decoupling scheme, we follow the stan¬ 
dard decoupling method^ and introduce three basic rules 
for the multi-orbital model as: (1) Some specific second- 
order Green’s functions, such as {{uiasCiaa \ 

{{nucrCiaa \ and {{nuaCiaa \ c]pcr)), are regarded 

as basic Green’s functions that receive the same treat¬ 
ment as the first-order Green’s function. That is, we also 
need to find the equations of motions of those Green’s 
functions. (2) In a third-order Green’s function with n 
operators, such as {{nu^'niaaCiaa \ c/pa))’ introduce 
the average value of {nucr') to get the lower-order Green’s 
function as {niicr'){{niag-Ciacr I because its orbital 

index I is different from a of the other operators. (3) We 
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do not consider the effects of the correlation functions 
that contain spin flip terms. As a result, the equations 


of motion of some basic second-order Green’s functions 
can also be obtained as 


(w -|- /ict 


(w -|- /Xq 


(w -|- /Xq 
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s^ol s^a. 
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b£Fe,m bGFe,m 
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+ Ui{nilg) {{riiaaCiaa \ ^jpa)) ~ ^ ^ t^fi{{nila) {{dba \ + 2{dlgCilg) {{Ciaa \ Cj^a)))' (B3) 

b^Cu 


* yunsong@bnu.edu.cn 

^ Y. Kamihara, H. Hiramatsu, M. Hirano, R. Kawamura, H. 
Yanagi, T. Kamiya, and H. Hosono, J. Am. Chem. Soc. 
128 , 10012 (2006). 

^ Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, J. 

Am. Chem. Soc. 130 , 3296 (2008). 

® X. Chen, P. Dai, D. Feng, T. Xiang, and F.-C. Zhang, Nat. 
Sci. Rev. 1, 371 (2014). 

C. W. Chu, Nat. Phys. 5, 787 (2009). 

® P J Hirschfeld, M M Korshunov, and 11 Mazin, Rep. Prog. 
Phys. 74 , 124508 (2011). 

® H. Alloul, J. Bobroff, M. Gabay, and P. J. Hirschfeld, Rev. 
Mod. Phys. 81, 45 (2009). 

^ A. V. Balatsky, I. Vekhter, and J.-X. Zhu, Rev. Mod. Phys. 
78 373 (2006). 

® N. Ni, A. Thaler, J. Q. Yan, A. Kracher, E. Colombier, S. 
L. Budko, P. C. Canfield, and S. T. Hannahs, Phys. Rev. 
B 82 , 024519 (2010). 

® J. Li, Y. F. Guo, S. B. Zhang, J. Yuan, Y. Tsujimoto, X. 
Wang, C. I. Sathish, Y. Sun, S. Yu, W. Yi, K. Yamaura, 
E. Takayama-Muromachiu, Y. Shirako, M. Akaogi, and H. 
Kontani, Phys. Rev. B 85, 214509 (2012). 

A. J. Williams, T. M. McQueen, V. Ksenofontov, C. Felser 
and R. J. Gava, J. Phys.: Gondens. Matter 21, 305701 
(2009). 

J. A. McLeod, A. Ruling, R. J. Green, T. D. Boyko, N. A. 
Skorikov, E. Z. Kurmaev, M. Neumann, L. D. Finkelstein, 
N. Ni, A. Thaler, S. L. Bud’ko, P. G. Canfield, and A. 
Moewes, J. Phys: Condens. Matter 24, 215501 (2012). 

Y. J. Yan, P. Cheng, J. J. Ying, X. G. Luo, F. Chen, H. 


Y. Zou, A. F. Wang, G. J. Ye, Z. J. Xiang, J. Q. Ma, and 
X. H. Chen, Phys. Rev. B 87, 075105 (2013). 

P. Cheng, B. Shen, F. Han, and H.-H. Wen, Eur. Phys. 
Lett. 104 37007 (2013). 

S. Ideta, T. Yoshida, 1. Nishi, A. Fujimori, Y. Kotani, K. 
Ono, Y. Nakashima, S. Yamaichi, T. Sasagawa, M. Naka- 
jima, K. Kihou, Y. Tomioka, C. H. Lee, A. lyo, H. Eisaki, 

T. Ito, S. Uchida, and R. Arita, Phys. Rev. Lett. 110 , 
107007 (2013). 

M. Merz, P. Schweiss, P. Nagel, Th. Wolf, H. v. Lohneysen, 
and S. Schuppler, arXiv: 1306.4222 

A. F. Wang, J. J. Lin, P. Cheng, G. J. Ye, F. Ghen, J. Q. 
Ma, X. F. Lu, B. Lei, X. G. Luo, and X. H. Chen, Phys. 
Rev. B 88, 094516 (2013). 

S. T. Cui, S. Kong, S. L. Ju, P. Wu, A. F. Wang, X. G. 
Luo, X. H. Ghen, G. B. Zhang, and Z. Sun, Phys. Rev. B 
88, 245112 (2013). 

T. -W. Huang, T.-K. Ghen, K.-W. Yeh, G.-T. Ke, C. L. 
Ghen, Y.-L. Huang, F.-G. Hsu, M.-K. Wu, P. M. Wu, M. 
Avdeev, and A. J. Studer, Phys. Rev. B 82 , 104502 (2010). 
T. Berlijn, C.-H. Lin, W. Garber, and W. Ku, Phys. Rev. 
Lett. 108 , 207003 (2012). 

T. Miyake, K. Nakamura, R. Arita, and M. Imada, J. Phys. 
Soc. Jpn. 79 , 044705 (2010). 

M. Nakajima, S. Ishida, T. Tanaka, K. Kihou, Y. Tomioka, 
T. Saito, G. H. Lee, H. Fukazawa, Y. Kohori, T. Kakeshita, 
A. lyo, T. Ito, H. Eisaki, and S. Uchida, Sci. Rep. 4 5873 
(2014). 

P. Werner, M. Casula, T. Miyake, F. Aryasetiawan, A. J. 
Millis, and S. Biermann, Nat. Phys. 8, 331 (2012). 





13 


D.-Y. Liu, Y.-M. Quan, D.-M. Chen, L.-J. Zou, and H. Q. 
Lin, Phys. Rev. B 84, 064435 (2011). 

L. de Medici, S. R. Hassan, M. Capone, and X. Dai, Phys. 
Rev. Lett. 102, 126401 (2009). 

L. de Medici, G. Giovannetti, and M. Capone. Phys. Rev. 
Lett. 112, 177001 (2014). 

Y. -M. Quan, L.-J. Zou, D.-Y. Liu, and H. Q. Lin, J. Phys. 
Conden. Matt. 24, 085603 (2012). 

K. Koepernik and H. Eschrig, Phys. Rev. B 59, 1743 
(1999). 

J. P. Perdew and Y. Wang, Phys. Rev. B 45, 13244 (1992). 
P. Soven, Phys. Rev. 156, 809 (1967). 

“ D. W. Taylor, Phys. Rev. 156, 1017 (1967). 

K. Koepernik, B. Velickj), R. Hayn, and H. Eschrig, Phys. 
Rev. B 55, 5717 (1997). 

M. Daghofer, A. Nicholson, and A. Moreo, Phys. Rev. B 
85 184515 (2012). 

M. Daghofer, A. Nicholson, A. Moreo, and E. Dagotto, 
Phys. Rev. B 81, 014511 (2010). 

T. A. Costi and A. Liebsch, Phys. Rev. Lett. 99, 236404 
(2007). 

K. Bouadim, G. G. Batrouni, and R. T. Scalettar, Phys. 
Rev. Lett. 102, 226402 (2009). 

H. Lee, Y.-Z. Zhang, H. O. Jeschke, R. Valenti, and H. 
Monien, Phys. Rev. Lett. 104, 026402 (2010). 

D. N. Zubarev, Soviet Phys. Usp. (English Trask) 3, 320 
(1960). 

J. Hubbard, Proc. Roy. Soc. A 276, 238 (1963). 

I. 1. Mazin, D. J. Singh, M. D. Johannes, and M. H. Du, 
Phys. Rev. Lett. 101, 057003 (2008). 

K. Byczuk, W. Hofstter, and D. Vollhard, Chapter 20, 50 
years of Anderson Localization ed. E. Abrahams (World 
Scientific, 2010). 

D. C. Licciardello and D. J. Thouless, Phys. Rev. Lett. 35, 
1475 (1975). 

C. He, Y. Zhang, B. P. Xie, X. F. Wang, L. X. Yang, B. 
Zhou, F. Chen, M. Arita, K. Shimada, H. Namatame, M. 
Taniguchi, X. H. Chen, J. P. Hu, and D. L. Feng, Phys. 
Rev. Lett. 105, 117002 (2010). 

L. Boeri, O. V. Dolgov, and A. A. Golubov, Phys. Rev. 
Lett. 101, 026403 (2008). 

K. Haule, J. H. Shim, and G. Kotliar, Phys. Rev. Lett. 
100, 226402 (2008). 

S. Zhou and Z. Wang, Phys. Rev. Lett. 105, 096401 (2010). 

M. Yi, D. H. Lu, R. G. Moore, K. Kihou, C.-H. Lee, A. 
lyo, H. Eisaki, T. Yoshida, A. Fujimori, and Z.-X. Shen, 
New J. Phys. 14 073019 (2012). 

Z. -R. Ye, Y. Zhang, B.-P. Xie, and D.-L. Feng, Chin. Phys. 
B 22, 087407 (2013). 


J. Knolle, V. B. Zabolotnyy, 1. Eremin, S. V. Borisenko, 

N. Qureshi, M. Braden, D. V. Evtushinsky, T. K. Kim, 

A. A. Kordyuk, S. Sykora, Ch. Hess, 1. V. Morozov, S. 
Wurmehl, R. Moessner, and B. Buchner, Phys. Rev. B 86, 
174519 (2012). 

L. X. Yang, Y. Zhang, H. W. Ou, J. F. Zhao, D. W. Shen, 

B. Zhou, J. Wei, F. Chen, M. Xu, C. He, Y. Chen, Z. D. 
Wang, X. F. Wang, T. Wu, G. Wu, X. H. Chen, M. Arita, 

K. Shimada, M. Taniguchi, Z. Y. Lu, T. Xiang, and D. L. 
Feng, Phys. Rev. Lett. 102, 107002 (2009). 

H. Ding, P. Richard, K. Nakayama, K. Sugawara, T. 
Arakane, Y. Sekiba, A. Takayama, S. Souma, T. Sato, T. 
Takahashi, Z. Wang, X. Dai, Z. Fang, G. F. Chen, J. L. 
Luo, and N. L. Wang, Eur. Phys. Lett. 83, 47001 (2008). 
G. Lee, H. S. Ji, Y. Kim, C. Kim, K. Haule, G. Kotliar, B. 
Lee, S. Khim, K. H. Kim, K. S. Kim, K.-S. Kim, and J. H. 
Shim, Phys. Rev. Lett. 109, 177001 (2012). 

J. Ferber, K. Foyevtsova, R. Valenti, and H. O. Jeschke 
Phys. Rev. B 85, 094505 (2012). 

L Y Xing, H Miao, X C Wang, J Ma, Q Q Liu, Z Deng, 
H Ding and C Q Jin, J. Phys.: Condens. Mat. 26, 435703 
(2014). 

S. F. Wu, P. Richard, A. van Roekeghem, S. M. Nie, H. 
Miao, N. Xu, T. Qian, B. Saparov, Z. Fang, S. Biermann, 
Athena S. Sefat, and H. Ding, Phys. Rev. B 91, 235109 
(2015). 

D. J. Singh, Phys. Rev. B 79, 153102 (2009). 

V. K. Anand, P. K. Perera, A. Pandey, R. J. Goetsch, A. 
Kreyssig, and D. C. Johnston, Phys. Rev. B 85, 214523 
( 2012 ). 

R. Kraus, V. Bisogni, L. Harnagea, S. Aswartham, S. 
Wurmehl, G. Levy, 1. S. Elfimov, B. Buchner, G. A. 
Sawatzky, and J. Geek, Phys. Rev. B 87, 134516 (2013). 

A. L. Efros and B. 1. Shklovskii, J. Phys. G: Solid State 
Phys. 8. L49 (1975) 

Y. Song, S. Bulut, R. Wortis, and W. A. Atkinson, J. 
Phys.: Condens. Matter 21, 385601 (2009). 

Q. Luo, G. Martins, D. X. Yao, M. Daghofer, R. Yu, A. 
Moreo, and E. Dagotto, Phys. Rev. B 82, 104508 (2010). 

B. R. Kubo, J. Phys. Soc. Jpn. 12, 570 (1957). 

P. A. Lee and D. S. Fisher, Phys. Rev. Lett. 47, 882 (1981). 
A. Wang, Single crystal growth and physical properties of 
iron-based superconductors with alkali metal, Ph.D. thesis, 
University of Science and Technology of Ghina, 2014. 

W. Nolting, Ghapter 4, Fundamentals of Many-body 
Physics: Principles and Methods, (Springer-Verlag, Berlin, 
2009). 



